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Abstract. In this paper it is shown that the characteristic vector field asso- 
ciated to a first order PDE: 

h(x, Vz(x)) = with x £ R n 

has the same form of an infinitesimal generator of an odd-symplectic transfor- 
mation with contact Hamiltonian h(x,p) on the level set h = 0. We also study 
under which condition such PDE has a characteristic vector field commuting 
with a generator of an odd-symplectic action. 



Introduction 

Very recently R. Cushman [0 presented the notion of odd-symplectic group. This 
notion was formulated in || but differently in Q ; more recently in j^, [| . Such group 
can be defined on vector spaces with odd dimension. Loosely speaking, the odd- 
symplectic group lies between the standard symplectic group and the contact group. 
More precisely, it interpolates symplectic groups in different dimensions according 
to the following exact sequence: 

. . . -> Sp(2n + 2, R) -» Sp{2n + 1, R) -> Sp(2n, !£)—►... 

Although the definition of the odd-symplectic group is natural, at the moment, 
there are few examples of applications. The aim of this paper is to show that 
odd-symplectic transformations rise naturally in the geometric theory of first order 
PDE and in particular in the study of the characteristics of PDEs of the following 
form: 

(1) h(x, Vz(x)) = 0, 

where x £ R™ h(., .) is a smooth function on R 2n . 

Here is an outline of the paper: in the first section the odd-symplectic group is 
defined. In the following sections we recall the theory of characteristics for first order 
PDEs and we show that the characteristic vector field associated to ([!]) coincides 
with an infinitesimal generator of an odd-symplectic transformation. We show also 
that 

h{x, Vz(x)) = - < gVz, Vz > + < v, Vz > +k 

with g constant n x n symmetric matrix, k £ R and v G R" has a symmetry 
whose infinitesimal generator is given by an element in the Lie algebra of Sp(2n + 
1,R). Finally we show that there is an isomorphism between certain class of first 
order PDE and the Lie algebra of the odd-symplectic group; moreover we show 
the connection of the odd-symplectic group with Hamilton-Jacobi theory and the 
Eikonal equation. 
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1. Odd-symplectic group and contact geometry 

Let us consider a symplectic vector space (V, to), dim V = 2n+2. The symplectic 
group is: 

(2) S P {V,uj) = {AeGl(V)\A*uj = uj}. 

The odd-symplectic group Sp(V,ui) Vo can be defined as follows: let vq £ V be a 
non-zero vector, then: 

(3) Sp{V,w) V0 ={AeSp(V,cj)\Av = v } 

Sometimes Sp(V t ui) Vo is denoted by Sp(2n + 1,R). This group has been studied 
in ||. It turns out that the odd symplectic group arises naturally in contact 
geometry. Let us follow M, Consider a contact manifold (A4, 9) (dim M. = 2n + l) 
where locally, in Darboux coordinates, the contact form is: 

(4) 6 = dz - yidx 1 

In |Q] it is shown that one can symplectify (M.,0) into {Af,oS) = (R+ x M,d(t9)) 
where t > is the coordinate on R + . Consider a one parameter group of symplec- 
tomorphisms F s : TV — > TV of the form: 

(5) F s (t,m) = (t,G s (t,m)). 

As shown in j| the infinitesimal generator of F s at t = 1, y = is an element of 
the algebra sp(2n + 1, R) of Sp(2n + 1, R). If A" is vector field on TV such that 

(6) L x (t) = and L x {d{t6)) = 
then its flow is a map of the form (|]) . Since 

L x (d(te)) = d(i x (d(te))) = o 

we deduce that (^) implies: 

(7) Lx(i) = and i x d(t6) = dh for some h : JV-> R 

To see what this means, let us use Darboux coordinates. If we denote 

dt d dz d dx l d dyi d 
_i_ _j_ _j_ zl 

h ds dt ds dz ds dx l ds dyi 



f =0 

1 dh 

ds t dyt 
dyj _ 1 dh 

ds i^dx z ^ 

dz dfi I y%_ dh 

ds dt t dyi 



then (0) gives: 
(8) 

with the condition: 

m f=» 

This has the following consequences: 

Remark 1.1. Since dt/ds — 0, the flow preserves t — const. Moreover 

L x% (h) = 0, 

the flow in M associated to the vector field preserves the (constant) Hamiltonian 
function h. 
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Remark 1.2. For t = 1 and h = t ■ h(x,y) one recovers the contact vector field 
defined in 0, that we call X\ which is equal to: 



(10) 



dx 1 
ds 

dyi 
ds 

dz 



dh 

dyi 

dh 

dx 1 - 

-h + 



Remark 1.3. The contactomorphisms generated by vector fields X\ form a sub- 
group of the group of contact transformation. In fact, in [^j V.Arnol'd defined a 
contact Hamiltonian vector field without requiring that the associated one parameter 
transformation has the form ^). Arnol'd's definition leads to a vector field that 
generates a contact transformation but does not preserve the contact Hamiltonian 
h. 



Remark 1.4. 

(11) 

1, 



Consider h u : M 



such that: 



< Vo, x > kt 

t 

where a, (3, 7 are nx n real matrices with (3 T = (3 and 7 T = 7. One verifies that the 
1 is the infinitesimal generator of the one parameter group 
expsu with u G sp{2n + 1,K), that is, 

exp(su) • (t, x, y, z). 



h u = ^[— < jx, x > +2t < ax, y > +t 2 < (3x, y >]+ < Vi,y > 



vector field X-^ at t 



X K 



ds 



(12) 
Here 
(13) 



where J is the standard symplectic structure in 
sp(2n, K) such that 




v <E 



k G K and A 6 



(14) 



JA = 



-7 

.T 



Putting (x,y) =w, a simple computation shows that on level set t = 1: 



(15) 



X 



t=l 



(v + Awf 



d_ 

dw 



(k + (Jv) T w) 



d_ 

Op. 



2. Theory of characteristics for first order PDE and 

ODD-SYMPLECTIC VECTOR FIELDS 



In this section we shall follow the presentation of first order PDE given in |l| . 
Let us consider the case in which the manifold Ai is the first jet bundle associated 
to R n . Locally such space can be modeled by a x G R n and a germ of a function 
in x considered up to its gradient. To each m G J 1 (E") we associate (x, z(x), V2) 
with z(.) a germ of a smooth function in K™. 
The manifold J 1 (i? ra ) is a contact manifold with contact 1-form 

9 = dz — yidx 1 

which vanishes on any germ z(.) such that dzjdx 1 — yi. The first jet space is the 
natural place to study the geometry of first order PDEs. In fact any PDE can be 
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understood as a submanifold of J 1 (i?") (see Let us suppose that we have a 
smooth function on J 1 (R n ): 

Bh 

h : J 1 (R n ) — ► R with — = 0. 

az 

h = defines a manifold in J 1 (i?™) which corresponds to the following PDE: 

(16) /i(a;,z(x), Vz(x)) = 

To equation ( |l6| ) we can associate a vector field X c whose integral curves satisfy 
(fl0|). We pose the natural question whether X c is related to the characteristic 
vector field, which converts a first order PDE into a system of ordinary differen- 
tial equations. In order to answer this question we need to recall the theory of 
characteristics. See Q for more details. 

2.1. Theory of characteristics for first order PDE. Let us assume that we 
have a first order PDE defined by the level set of h : J 1 (i?™) — » R. In the first jet 
space J 1 (R") we can introduce contact plane TL X at any x E R™ spanned by the 
vectors V 

(17) i v = 

Equation ( |l7| ) determines a codimension 1 submanifold E 2n . Its tangent space 
TE 2n can be described by the vectors V such that: 

(18) i v dh = 

Definition 2.1. T/ie submanifold E 2n is called non- characteristic, if TE 2n is 
transversal to the contact plane If. 

Definition 2.2. The intersection of vector spaces 

p x = T x E 2n n n x 

is called characteristic plane at x. Therefore P is defined by all vectors V such that 



(19) i v 9 = and i v dh = 0. 

We shall consider only the case in which E 2n is not characteristic. Note that 
since TE 2n and II lie in J 1 (R n ) we have: 

dimP^ = dimTEl 71 + dimll x - dimTJ 1 (R n ) = In- 1 

The 2-form dd, associated to the contact form 8, defines a skew product in TJ 1 (R n ). 
This allows us to select a vector in each characteristic plane P x which is skew- 
orthogonal to P. The characteristic vector field for ( |l6| ) is given by X c : 

(1) ix c — 0, X c is a contact vector field, 

(2) d9(X c , V) = for all V E P. 

In symplectic geometry one can show that the skew-orthogonal complement to a k 
dimensional vector space is a 2n — k vector space. The characteristic plane P has 
dimension In — 1. Therefore the characteristic direction X c is uniquely determined. 
For a given h : J 1 (R n ) — * R, in local coordinates the vector field X c has components: 



(20) 



dx l dh 

ds dyi 
dyj _ dh 

dz on 

ds y % dyi 



y* d dz 
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2.2. Cauchy Problem. The theory of characteristics is used to solve the Cauchy 
problem. For an equation 

h(x, z(x), Vz(i)) = 

Cauchy data are: a manifold 7 C K" and the values z\ — 4>{x). We have seen 
that the equation h = corresponds to submanifold E 2n in J 1 (R"). Cauchy data 
define a submanifold N C E n called the manifold of initial data. It turns out that 
the Cauchy problem is solvable (locally) if N is not characteristic, that is if the 
projection of the characteristic direction to R™ is transversal to the tangent plane 
to N (see [Q). Given a parametrization x 1 {\) of 7, it is possible to construct the 
family of initial data: 

x = x 1 (X),z = (j>(x y (X)),po =PoW 
where Po(X) is such that: 

h(x^X),4>{x^(X)),p (X)) =0 

which forms the family of initial data for the characteristic vector field (|20|). 

Remark 2.1. The characteristic vector field ( pZ^ does not depend on whether the 
PDE is given by h = or h — c with c ^ 0, but the Cauchy problem does. In fact 
the value of the level set of h is contained in pq. 



3. The main result 

We can state: 

Theorem 3.1. For all the PDEs represented in J (R n ) by a smooth function h : 
J l {R n ) — » R. of the form: 

h{x,Vz{x)) = 0, 



the odd-symplectic Hamiltonian vector field commutes with the characteristic 
vector field (\2w. Moreover on h — 0, the two vector fields coincide. 



Proof. Consider the function h : J^R") -> R as a contact Hamiltonian. Then by 
inspection one finds that the infinitesimal generator can be written as follows: 

(21) Xl=X c -h^-. 

Since h dos not depend on z, therefore X c (h) = 0. Also the Lie bracket: 

d 



[X c ,X h ] 



X c ,h— 
oz 



vanishes. 

Now we want to compare X\ (JToj) and X c @. We notice that if h : J 1 (R n ) -> R 
does not depend on z, the first two components are equal. The component dz/ds 
differs because of the term — h in (|l0|). We have 

L x i(h) = 

that is h is preserved by the flow of X\ . Therefore in the manifold h = the vector 
field X^ is equal to the characteristic vector field X c □ 
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3.1. First order PDE and sp(2n + 1,R). Let us consider a particular class of 
first order PDEs in R™, namely, 

dz(x) dz(x) , dz(x) .- t j dz(x) , , 

or more compactly: 

(23) < BVz(x), Vz(x) > + < e, Vz(ir) > + < /, x >= h 

Then we have: 

Proposition 3.1. For a given first order PDE of the form (1131), if 



1 f 2a b 

2 { b T 2c 



(where . T is the transpose) is symmetric, then we can associate u £ sp(2n + 1,R) 



(24) 

where k = ho, Jv — (e, /) and —JB = A, to the characteristic vector field of f 
If a = b = and v = (i>i,0) in B then the characteristic vector field X c associated 
to commutes u. 




Proof. The proof is essentially based on theorem 3.1 and on remark 1.4. According 
to theorem 3.1, equation ( p3| ) can be interpreted as a contact Hamiltonian h(x,y) 
where y — Vz(i), and the contact vector field on h = is equal to the characteristic 
vector field. If the conditions on the matrices a, 6, c are satisfied the PDE (p3|), on 
J 1 (R") gives rise to a contact Hamiltonian of the form: 

h(x, y, 1) =< w, J Aw > + < Jv, w > +ho, where w — (x, y) 



According to remark 1.4, we can associate to such Hamiltonian an element of 
sp(2n + 1,R) of the form (§J). 
Let us consider the Hamiltonian 

(25) h u —< ax, x > +2t < bx, y > +t 2 < cx, y > + < v\,y > — < v%, x > — — hot 
The corresponding infinitesimal generator X~ is: 

where 

(27) =<vi,y > - < v 2 , x > +h 
The characteristic vector field X c of ( p3| ) is 

(28) X c = ti + Z c -?-, 

oz 

where 

(29) i = 2{bx + cy)-^--2{ax + by)-^- 

(30) Z c = 2 < bx,y > +2 < cy,y > + < v x ,y > 
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Since h does not depend on z, the commutator is given by: 

[xi h) x c } = (az c )-m))^- z 

After some algebra one finds: 

£(Z C ) - £(Z$) = 4<bx, b T y > +2 < cy, b T y > + < bx, v 2 > + < cy, v 2 > + 

(31) -2 < bx, ax > -2 < bx, by > -4 < ax, cy > -4 < by, cy > 

Last expression is identically zero for every (x, y) only for a = b = and t> 2 = 0. □ 

3.2. Hamilton-Jacobi theory. In certain cases in Hamiltonian mechanics the 
construction of coordinate transformations that allows to integrate explicitly the 
equations of motion is based on the solution of Hamilton-Jacobi equation. We refer 
the reader to Q. Suppose that we have a Hamiltonian system on (R 2n ,w). The 
Cartan one-form, written in local coordinates, is: 

(32) 6 =Pidq l - H{p,q,t)dt 

where H : R 2 " x M — > R in the Hamiltonian. is an integral invariant for the 
Hamiltonian flow. Consider the extended phase space A4 = M. 2n x R and two set 
of coordinates: 

(q,p,t) and (x,y,r) 
There are two functions K(x,y,r) and S such that 

(33) Pidq 1 — H(p, q, t)dt = yidx 1 — K(x, y, r)dr + dS 

The function S is called a "generating function" because it allows one to define a 
canonical transformation. 

Equation (|3^) can be used to find the canonical transformation which realizes 
(q,p,t) — » (x,y,r). This transformation is usually consider in the case t = r. 
In this case, by means of Legendre transformation one can write (k3^ as: 



(34) pidq 1 - H(p, q, t)dt = x l dy { - K(x, y, t)dt + dS{q, y, t). 

Then one usually chooses K — and obtains: 



(35) 



w = v . 

dq* ~ ft 



tf(>!f) + f = 



The third equation of ( |35| ) is the Hamilton-Jacobi equation and determines the func- 
tion 5* which completely describes the canonical transformation (q,p,t) — > (x,y,t). 
The Hamilton-Jacobi equation is defined in J 1 (M n+1 ) and fulfills the conditions of 



theorem 3.1 



3.3. Eikonal equation. In geometric optics and WKB approximation of quantum 
mechanics WKB approximation the Eikonal equation is relevant. This is a first 
order PDE given by: 

(36) h i ^^-N(x)=0 
y ' 2 dx l dxi v ' 

S represents the wave front and N is a function: 

N(x) = (n(x)/c) 2 , 
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G = 1 ° 



where in optics is square ratio of the refraction index n(x) to the speed of light c; 
where in WKB theory: 

N(x) = V(x) - E, 
namely, potential minus the in energy. 

When in optics we consider a medium with a piece-wise constant refraction index 
(multilayers), N(x) can be easily approximated by piece- wise constant function. 
In quantum mechanics this happens when the potential is piece- wise constant. In 
these cases proposition [O] can be applied in such layer. For instance if 

N(x) — n^/c 2 for x in the A-layer 

then in fc t?l -layer we have symmetry generated by 

(37) 

with 

<^ — 

2 V 1 

A similar analysis can be carried out in the case of semi-classical description of 
system which have piece-wise constant potential e.g. multi-well potential. 
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